Enhancing the Detection of Natural Thermal Entanglement with Disorder 
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Physical systems have some degree of disorder present in them. We discuss how to treat natural, 
thermal entanglement in any random macroscopic system from which a thermodynamic witness 
bounded by a constant can be found. We propose that functional many-body perturbation theory 
be applied to allow either a quenched or an annealed average over the disorder to be taken. We find 
when considering the example of an XX Heisenberg spin chain with a random coupling strength, 
that the region of natural entanglement detected by both witnesses can be enhanced by the disorder. 

PACS numbers: 03.67.-a 05.40.-a 61.43.-j 



Many-body systems have been of great interest to the 
condensed matter community for numerous years. More 
recently, ideas from quantum information have been ap- 
pUcd to such systems, allowing new methods of quan- 
tifying their properties to be utilised. Entanglement 
in many-body systems, [ij is of particular significance 
due both to fundamental interest and applications in 
quantum computation. Approaching the thermodynamic 
limit and at finite temperatures, this natural entangle- 
ment becomes relevant in real, macroscopic systems ^. 

Although entanglement between any two qubits, or 
even subsystems, in a large system can be (compara- 
tively) easily calculated, quantifying entanglement be- 
tween more qubits is much more complex [l|. However, 
using an entanglement witness we can still find critical 
values of the system's parameters below which the system 
is entangled. An entanglement witness is an observable 
whose expectation value is bounded for any separable 
state. Using a witness is advantageous as we can de- 
tect natural, and therefore thermal, entanglement of all 
classes while other measures detect only bipartite entan- 
glement and/or can be used only at zero temperature. 

When an entanglement witness can be derived from 
thermodynamic quantities, in particular the partition 
function, we refer to it as a thermodynamic witness [3|. 
If a system contains some randomness, an average over 
this disorder must be performed. The method used to 
calculate the average is dependent on the time taken to 
measure the witness in an experimental situation, giving 
either a quenched or an annealed average [^. For each 
type of average, the witness will detect a different region 
of entanglement. It is therefore important to distinguish 
which we are able to measure experimentally. Previous 
research into average entanglement [^, m 01 has concen- 
trated on spin chains and used sampling, which corre- 
sponds to calculating a quenched average numerically, or 
a renormalisation group approach to calculate how the 
randomness in the chain affects entanglement. These 
methods cannot take the differences between quenched 
and annealed entanglement into account. In addition. 



these papers look at zero temperature, bipartite entan- 
glement. It is both in these areas, and by introducing a 
perturbative technique to deal with disorder in the con- 
text of entanglement, that our work offers new insight. 

This letter studies natural entanglement in random 
systems. After averaging over the disorder, the regions 
of entanglement detected by the quenched and annealed 
witnesses can be related by Jensen's inequality [8|. As 
these averages are difficult to perform on a large system, 
we introduce the concept of using the well established 
many-body perturbation theory to calculate the thermo- 
dynamic entanglement witness. This method is very gen- 
eral and can be used, at least numerically, for any system. 

To exemplify how these different averages affect a sys- 
tem's entanglement, we consider an XX spin half chain 
in a uniform magnetic field with an extra random cou- 
pling strength. This system in particular is used as 
the majority of the calculation can be done analytically 
rather than numerically. Hence we calculate both the 
quenched and annealed witnesses for this system using 
many-body perturbation theory and plot the results. We 
find that the randomness can both enhance and create 
regions of entanglement detected by both witnesses. We 
stress that it is the region of entanglement detected by 
the witness that is enhanced by the disorder, rather than, 
necessarily, the entanglement itself. Intuitively however, 
increasing the disorder can increase the entanglement it- 
self, hence wc discuss possible mechanisms for this. 

In order to calculate an entanglement witness in a ran- 
dom system, we must perform an average over the disor- 
der. As discussed above, the appropriate average to use 
is dependent on how fast the measurement process is in 
comparison to the speed of the fluctuations of the vari- 
able in the system. As different experimental conditions 
require different averages to be taken of the system, it is 
important to determine which is used. 

First we consider the possibility that we can measure 
each state of the system faster than the disorder fluctu- 
ates. This is called a quenched average. To calculate this, 
we first perform a thermal average, and then average over 



the random variable. Rather than calculating the proper- 
ties of the system, we calculate an average of all possible 
systems. Otherwise, if the system fluctuates faster than 
we can measure it, we cannot measure each state the sys- 
tem takes, and we should use an annealed average. Here, 
the system thermalises before we can measure the state; 
hence the average over the disorder and the thermal av- 
erage are taken at the same time. Experimental studies 
on high temperature superconducting materials discuss 
situations when each of these averages occur [91 . 

As each type of average relates to a different standard 
of measurement, each will give a different region of en- 
tanglement. It is therefore important to clarify which 
is needed when measuring a system. The regions of en- 
tanglement detected by the quenched and annealed av- 
erages can be related by invoking Jensen's inequality, 
{f {X)) < f{X), where is / is a concave function act- 
ing on X and the average, (• • •) is over the disorder. For 
a convex function, the inequality is reversed. The left 
hand side of the inequality is the quenched average, and 
the right hand side is the annealed average. We are in- 
terested in how these averages can change the amount of 
entanglement we can measure, hence we must consider 
how wc can use entanglement witnesses in this context. 

For any system of finite dimension for which we can 
calculate a thermodynamic witness with some constant 
value, a, as the bound for detecting entanglement, we 
can apply the following argument. If the witness, W, is 
a convex function, the entangled region detected which 
we define as E = a — W is always concave. Similarly, E 
is convex if the witness is concave. Hence the application 
of Jensen's inequality to cither scenario tells us which 
average allows us to detect the most entanglement. 

Although entanglement has previously been studied 
using thermodynamic witnesses calculated from the par- 
tition function of the system, to discover how disorder 
affects thermal, macroscopic (natural) entanglement, we 
instead use functional many-body perturbation theory. 

We set out an outline of the method for a fermionic sys- 
tem below. Due to the Jordan- Wigner transformation, 
this also applies to any spin half system. Bosonic sys- 
tems can also be studied using this formalism, but com- 
plex numbers are used in place of Grassmann variables. 
We first define a generating functional. Z[fj,r]]\fj^ri=o = Z 
where Z = Tr{e^^^) is the partition function of the sys- 
tem and ?7 is a Grassmann variable. Z[fj,i]] allows us to 
calculate expectation values using functional derivatives 
from an appropriate generating function, gf. 
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(1) 



where /? is the inverse temperature, H = Hq + Hi{vi) 
is the Hamiltonian of the system and vi is a random 
variable, with I labelling position. T is the time ordering 
operator. After putting ([T]) into path integral representa- 
tion [lO[ , we can take an average over the random variable 



of Z[?7, 77]". For clarity, we consider the random variables 
to be taken from a Gaussian distribution. 
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(2) 
where A is the variance. After calculating this average, 
we perform transformations to diagonalise Hq- The av- 
erage can then be rewritten as 
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where Zq is the partition function of the unperturbed 
system, A^^'^jjr, j-] is the perturbed term, and acts on 
Sf[Vi vli the generator re- written in terms of a Greens 
function. Expanding e >"*, we can use Feynman di- 
agrams to calculate {Z"'[fj,T]]). Functional derivatives 
can then be used to calculate an entanglement wit- 
ness for both the quenched and the annealed averages. 
We can then implement the replica trick, (InZ) = 
lim„^o -((^") ~ 1) in order to calculate the quenched 
witness. The annealed average is found simply by set- 
ting 71 = 1. The replica trick corresponds to considering 
n identical, non-interacting replicas of the system. Per- 
forming an average over the disorder of these replicas 
allows us to use A, the variance rather than a differ- 
ent random variable at each /. Rather than considering 
the experimental situation of measuring the system many 
times and performing an average, the replica trick allows 
us to instead use many copies of the same system. 

Perturbation theory is valid when the matrix elements, 
{m\Hi{vi)\n) , of the interaction term, Hi{vi), where \n) 
and \m) {rn ^ n) are possible eigenvectors of TJq, are 
much smaller than the energy difference between n and 
771. To ensure this condition is satisfied, we require t;; <C 1 
with high probability, and hence A ^ 1. 

To make the arguments above explicit, we calculate the 
average witness for a specific, predominantly analytically 
solvable, example: an XX spin half chain with a random 
coupling strength. The unperturbed Hamiltonian is 
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where J is the coupling strength between nearest neigh- 
bours, and i? is a uniform magnetic field. Considering 
the thermodynamic limit, iV — s- cxd, the spin chain can 
be diagonalised |ll| using a Jordan- Wigner transforma- 
tion ai = Hm^i '^m ® i'^f + ^'^f) /2 and a Fourier trans- 
form, a, = /:^ ||e^«'d(g) giving Ho = j% ^(Jcosq - 
B)S{q)d{q). The partition function is then 
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FIG. 1: This plot shows the values of uniform magnetic field, 
B and temperature, T, with J = 1 below which entanglement 
is detected by the quenched witness for various values of the 
disorder, A. Details of this plot are discussed in the text. 



FIG. 2; This plot shows the values of uniform magnetic field, 
B and temperature, T, with J = 1 below which entanglement 
is detected by the annealed witness for various values of the 
disorder, A. Details of this plot are discussed in the text. 



We now introduce the Hamiltonian H = Hq + Hi{iji) and 
consider Hi{ji) = -\ YaJ^^ ji {crfcrf+i + c^fcrf+i) where 
ji is a site dependent random coupling strength between 
nearest neighbours. Each random variable is taken from 
a Gaussian distribution centred at zero with variance A. 
To detect entanglement, we use an entanglement wit- 



ness 



12| which is proportional to the sum of the x and 



y components of the susceptibility, and is therefore mea- 
surable in practice [13|. 
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The system is entangled when the witness is greater than 
1. This bound is found using the Cauchy-Schwarz in- 
equality and properties of the density matrix IJ, [l5| . 

The entangled region, E = a ~ W, is concave as W is 
always convex due to the absolute sign. Using Jensen's 
inequality, this means the region of entanglement de- 
tected by the quenched witness is equal to or less than 
the region detected by the annealed witness. 

In order to use functional many-body perturbation 
theory, we must first define an appropriate generating 
function with which to calculate the witness. Letting 
r/3 



9f[v,v] = /o drj^i [m{T)ai{T)+T]i{T)al{i 
diagonalising the witness, 
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1 for the 



Srib{T',q') 

with a; = for the quenched average and x 
annealed average. 

Considering Hi{ji), we average over ji in ^ and 
calculate equation ^ from In Zq as given in ^ , AS'JJjj = 
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/o c^T-rfr' /^^ g Y.ab Vair, q)Gab{T - t', q)iib{T', q). 

After expanding e^'^'"* and calculating Feynman dia- 
grams to first order, we find the witness is 



W, 



TT 



dq cos q 



l{q) - AGW 



(8) 



with ave indicating either a quenched, q 
or annealed, a average. From the method 

we 



outlined above, 
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ieip)-e{q)f 

^n{q){l - n{q)) COS q J^^ dp n{p) COS p. 
^Peix) j^ ly'^ g^^^ g^^;) = 2Jcosa; - 25. 

As discussed previously, to ensure our use of perturba- 
tion theory is valid, we restrict the value of the disorder 
to A < 0.0001. We note that all the features found in 
Figs. H]) and ^ are present even when A <C 0.0001. 

The entanglement detected by the witness in an unper- 
turbed XX spin half chain, plotted in [l^l (Fig. (1) at 
6 = 0), exists below a critical uniform magnetic field, B, 
and temperature, T, and is below the coupling strength 
between nearest neighbour sites, J. Much attention has 
been focused on the study of quantum phase transitions 
[l6| and we note that Hq has one at B = J. The spec- 
trum of the unperturbed spin chain [17| at zero temper- 
ature is as follows: at S = 0, the ground state of the 
spin chain is a fully symmetric state of all permutations 
of an equal number of up and down spins. On increasing 
B, the number of up spins increases until, when B > J, 
they are all spin up and the system is in a product state. 



Fig. Il]) shows how the entanglement detected by the 
quenched witness behaves. Introducing the random cou- 
phng strength enhances the region detected by the wit- 
ness. Increasing A increases the entangled region de- 
tected compared to that in the original XX model. A 
higher variance increases the probability of a high value of 
ji (though this may be positive or negative) being picked 
out of the Gaussian distribution. As we consider small 
A and a larger coupling strength increases entanglement, 
so does increasing A. Further, a region of entanglement 
close to i? = J is created at low temperatures which 
again increases with A. We note that this extra region 
of entanglement is created at the quantum phase transi- 
tion of i^o- The annealed witness is plotted in Fig. ([2]) 
and gives a larger entangled region than the quenched 
witness everywhere except at T = where they coin- 
cide. Hence the plots demonstrate the importance of an 
accurate determination of how we measure the system 
experimentally. The behaviour of the plot as A changes 
remains the same as for the quenched witness. 

Alternatively, we can view the system in the following 
way: even if the ground state of a system is separable, 
there is a possibility of creating entanglement by increas- 
ing the temperature [18|. A finite temperature allows 
higher energy levels to mix with the ground state. If 
these newly available energy levels are themselves entan- 
gled, on mixing with a separable ground state, entangle- 
ment can be created and/or increase. Similarly, intro- 
ducing disorder into a system allows higher energy levels 
to be populated even at zero temperature, allowing for 
the possibility of creating or enhancing entanglement. 

We have also considered the effect of a random mag- 
netic field, bi on the witness. The plots for the quenched 
and the annealed witness have similar features as for jj. 
The region of entanglement close to B = J increases 
with A similarly to the random coupling strength plots. 
At lower B, increasing A decreases the entangled re- 
gion, contrary to the random coupling strength case. A 
random magnetic field increases the chance of spin flip 
events. For lower B this can only act to decrease the en- 
tanglement due to the highly entangled structure of the 
ground state of Hq. At higher B, particularly a,t B = J, 
this same effect will increase the entanglement due to the 
almost separable ground state at this point. 

This letter has studied entanglement in the limit of 
small disorder. The opposite, large disorder limit is also 
of interest; for bi, it is clear that large disorder would 
cause the system to become separable as H ^ '^^i^i ■ 
However, the case for ji is much more complex. Using 
the path integral formalism, we can consider Hi{ji) to 
be the unperturbed part of the Hamiltonian, and Hq to 
be the perturbation. As this cannot be done analytically, 
we do not pursue this here. We do, however, still expect 
that in the limit of large disorder the system should be- 
come completely mixed, and hence separable, though the 
effect of averaging over this disorder on the behaviour of 



the witness is nontrivial. This uncertainty is due to the 
equally high probability of a large positive, and hence 
entanglement enhancing, or negative, and hence entan- 
glement destroying, value of j;. While beyond the scope 
of this letter, this would be an interesting topic to study. 

We have discussed how the quenched and annealed 
entanglement witnesses are related both experimentally 
and via Jensen's inequality, and found that we can apply 
functional many-body perturbation theory to a random 
system in order to calculate these witnesses. We have 
found that randomness in a system can enhance the en- 
tanglement detected by a witness. In particular, we have 
considered a spin chain with a random coupling strength. 
This result is important in the experimental realisation of 
quantum systems due to entanglement being a resource 
in quantum computation as it demonstrates that disorder 
may actually be beneficial even in thermal systems. 
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